Scaling of Entanglement close to a Quantum Phase Transitions 
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In this Letter we discuss the entanglement near a quantum phase transition by analyzing the 
properties of the concurrence for a class of exactly solvable models in one dimension. We find that 
entanglement can be classified in the framework of scaling theory. Further, we reveal a profound 
difference between classical correlations and the non-local quantum correlation, entanglement: the 
correlation length diverges at the phase transition, whereas entanglement in general remains short 
ranged. 



Classical phase transitions occur when a physical sys- 
tem reaches a state below a critical temperature char- 
acterized by a macroscopic order [jlj. Quantum phase 
transitions occur at absolute zero; they are induced by 
the change of an external parameter or coupling con- 
stant [Q, and are driven by fluctuations. Examples in- 
clude transitions in quantum Hall systems 0, localiza- 
tion in Si-MOSFETs (metal oxide silicon field-effect tran- 
sistors; Ref. B) and the superconductor- insulator tran- 
sition in two-dimensional systems |5]](|. Both classical 
and quantum critical points are governed by a diverg- 
ing correlation length, although quantum systems possess 
additional correlations that do not have a classical coun- 
terpart. This phenomenon, known as entanglement [Q, 
is the resource that enables quantum computation and 
communication @]. The role of the entanglement at a 
phase transition is not captured by statistical mechan- 
ics - a complete classification of the critical many-body 
state requires the introduction of concepts from quan- 
tum information theory 0. Here we connect the the- 
ory of critical phenomena with quantum information by 
exploring the entangling resources of a system close to 
its quantum critical point. We demonstrate, for a class 
of one-dimensional magnetic systems, that entanglement 
shows scaling behaviour in the vicinity of the transition 
point. 

There are various questions that emerge in the study of 
this problem. Since the ground state wave-function un- 
dergoes qualitative changes at a quantum phase transi- 
tion, it is important to understand how its genuine quan- 
tum aspects evolve throughout the transition. Will en- 
tanglement between distant subsystems be extended over 
macroscopic regions, as correlations are? Will it carry 
distinct features of the transition itself and show scaling 
behaviour? Answering these questions is important for a 
deeper understanding of quantum phase transitions and 
also from the perspective of quantum information theory. 
So results that bridge these two areas of research are of 
great relevance. We study a set of localized spins coupled 
through exchange interaction and subject to an external 
magnetic field (we consider only spin- 1/2 particles), a 



model central both to condensed matter and information 
theory and subject to intense study fic[ |. 
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FIG. 1. The change in the ground state wave-function in 
the critical region is analyzed considering d\C(l) as a func- 
tion of the reduced coupling strength A. The curves corre- 
spond to different lattice sizes N = 11, 41, 101, 251, 401, oo. 
We choose N odd to avoid the subtleties connected with 
boundary terms jl6| . On increasing the system size, the mini- 
mum gets more pronounced, also the position of the minimum 
changes and tends as N~ 1,87 (left inset) towards the critical 
point A c = 1 where for an infinite system a logarithmic di- 
vergence is present (see equation (^|)). The right inset shows 
the behaviour of the concurrence C(l) itself for an infinite 
system. The maximum that occurs below A c is not related 
to the critical properties of the Ising model. As explained 
in the text, it is the change in the ground state and not the 
wavefunction itself that is a good indicator of the transition. 
The structure of the reduced density matrix, necessary to cal- 
culate the concurrence follows from the symmetry properties 
of the Hamiltonian. Reality and parity conservation of H to- 
gether with translational invariance fix the structure of p to 
be real symmetric with pn, p 2 2 = P33, P23, P14, pa as the 
only non-zero entries. 

O'Connors and Wootters showed that in the Heisen- 
berg chain the maximization of the entanglement at zero 
temperature is related to the energy minimization [ pr| . 
It is known that Werner states [O] can be generated 
in a one dimensional XY model [|13[ and that temper- 
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ature and magnetic field can increase the entanglement 
of the systems as shown for the Ising and Heisenberg 
models 0,|l5| ■ Finally we mention the study on the role 
of the entanglement in the density matrix renormaliza- 
tion group flow and the introduction of entanglement- 
preserving renormalization schemes |l6| ] . Here we address 
the problem of the relation between macroscopic order, 
classical correlations, and quantum correlations. There- 
fore we analyze the entanglement near the critical point 
of the XY model in a transverse field. Because of the uni- 
versality principle - the critical behaviour depends only 
on the dimension of the system and the symmetry of the 
order parameter - our results have much broader validity. 
We find that in the vicinity of a quantum phase transition 
the entanglement obeys scaling behaviour. On the other 
hand this analysis provides a clear distinction between 
the role of entanglement and correlations in quantum 
systems close to a critical point. (We have been made 
aware that similar work to that reported here is being 
performed by T. Osborne and M. Nielsen; T. Osborne, 
personal communication. p3) 



dimension- less coupling constant A = J/2h. For 7 = 1 
Eq. (Q) reduces to the Ising model whereas for 7 = it is 
the XY model. For all the interval < 7 < 1 the models 
belong to the Ising universality class and for N = 00 they 
undergo a quantum phase transition at the critical value 
A c = 1. The magnetization (a x ) is different from zero for 
A > 1 and it vanishes at the transition. On the contrary 
the magnetization along the z-direction (a z ) is different 
from zero for any value of A. At the phase transition the 
correlation length £ diverges as £ ~ |A — \ c \ v with v = 1 
(Refs. jl8| and p9|). We confine our interest to the en- 
tanglement between two spins, on position i and j, in 
the chain. All the information needed is contained in the 
reduced density matrix p(i,j) obtained from the ground 
state wave-function after all the spins except those at 
positions i and j have been traced out. The resulting 
p(i,j) represents a mixed state of a bipartite system; a 
good deal of work has been devoted to quantify the en- 
tanglement in this case 0,^0). We use the concurrence 
between sites i and j, related to the "entanglement of 
formation" , and defined as JItJ 
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FIG. 2. The finite size scaling is performed for the 
case of logarithmic divergences |2l)| . The concurrence, con- 
sidered as a function of the system size and the coupling con- 
stant, is a function of N 1/, "(\ — A m ) only, and in the case of 
log divergence it behaves as d\C(l)(N, A) - d\C(l)(N, Ao) ~ 
Q[N 1 ^(X-\ m )] - Q[N 1/l, (\o-X m )] where A is a non-critical 
value and Q(x) ~ Q(oo)lna; (for large x). All the data from 
N — 41 up to N = 2701 collapse on a single curve. The crit- 
ical exponent is v = 1, as expected for the Ising model. The 
inset shows the divergence of the value at the minimum as 
the system size increases. 

The system under consideration is a spin-1/2 ferromag- 
netic chain with an exchange coupling J in a transverse 
magnetic field of strength h. The Hamiltonian is 

j N N 

H = — 5>-y)of *£n + (i + 7)^i+i-fc£°i (1) 

i=l i=l 

where cr° are the Pauli matrices (a = x, y, z) and N is 
the number of sites. We assume periodic boundary con- 
ditions. It is convenient, for later purposes, to define a 



C(i,j) = max{0, riii, j)-r 2 (i, j)-r 3 (i, j)-r 4 (i, j)} (2) 

In Eq.(D) r a (i,j) are the square roots of the eigen- 
values of the product matrix R — p(i,j)p(i,j) in de- 
scending order; the spin nipped matrix p is defined as 
p = <r v (g) a v p*a y <E) cr y . In the previous definition, the 
eigenstates of <r z E {| f) | J.)} should be used. Trans- 
lation invariance implies that C(i,j) = C(\i — The 
concurrence will be evaluated as a function of the relative 
position between the spins and the distance |A — A c | 
from the critical point. The structure of the reduced den- 
sity matrix is obtained by exploiting symmetries of the 
model (see caption of Fig. 1). The non-zero entries of p 
can then be related to the various correlation functions, 
and the concurrence of the ground state is evaluated ex- 
actly starting from the results in Refs. |Q, p9| , and pi) . 

First we look at the Ising Model (7 = 1). The first 
question we consider is the range of the entanglement 
£e, that is, the maximum distance between two spins at 
which the concurrence is different from zero. The result is 
surprising: even at the critical point, where spin-spin cor- 
relations extend over a long range (the correlation length 
is diverging for an infinite system), the concurrence van- 
ishes unless the two sites are at most next-nearest neigh- 
bors. The truly non-local quantum part of the two-point 
correlations is nonetheless very short-ranged. 

In order to quantify the change of the many-body 
wave- function when the system crosses the critical point, 
we look at the derivatives of the concurrence as a function 
of A. In this case we need to consider only the nearest 
neighbour and next nearest neighbour concurrence. We 
first discuss the behaviour of the nearest neighbor con- 
currence. 



2 



0.014 ^-T 




FIG. 3. As in the case of the nearest neigh- 

bour concurrence, data collapse is also obtained for the 
next-nearest-neighbour concurrence C(2). In the figure, data 
for system size from N = 41 to N = 401 are plotted. The 
inset shows a peculiarity of the Ising model: C(2) has its 
maximum precisely at the critical point for arbitrary system 
size (note that the maximum decreases as the system size in- 
creases). Therefore we consider the second derivative to per- 
form the scaling analysis. It can also be seen that C(2) is two 
orders of magnitude smaller than C(l). For the smallest sys- 
tem sizes the concurrence is different from zero for \i — j\ —3 
and A > 1.05 (for N = 7; for N > 9, C(3) = for all A). 
In contrast the correlation functions are long-ranged at the 
critical point. 

The results for systems of different size (including the 
thermodynamic limit) are presented in Fig.^. For the in- 
finite chain d\C(l) diverges on approaching the critical 
value as 

d x C(l) = In IA-AJ + const. (3) 

Equation (|^) quantifies non-local correlations in the crit- 
ical region. One aspect of this system, particularly rel- 
evant for quantum information is the study of the pre- 
cursors of the critical behaviour in finite samples. This 
study is known as finite size scaling Q. In Fig.g the 
derivative of (7(1) respect to A is considered for different 
system sizes. As expected, there is no divergence for fi- 
nite N, but there are clear anomalies. The position of the 
minimum A m scales as A m ~ A c + Af -1 ' 87 and its value 
diverges logarithmically with increasing system size as 

d x C(l)\ Xm = -0.2702 In AT + const. (4) 

According to the scaling ansatz p2| in the case of log- 
arithmic singularities, the ratio between the two prefac- 
tors of the logarithm in Eqs.(|3|) and (Q) is the exponent 
which governs the divergence of the correlation length £ 
|||. In this case (8/3tt 2 « 0.2702) it follows that v = 1, 
as it is known from the solution of the Ising model Jl8| . 
By proper scaling ]22| ] and taking into account the dis- 
tance of the minimum of (7(1) from the critical point, it 



is possible to make all the data for different N collapse 
onto a single curve (Fig.||). 
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FIG. 4. The universality hypothesis for the entanglement 
is checked by considering the model Hamiltonian, defined in 
Eq.ffl), for a different value of 7. In this case we chose 7 = 0.5 
and N ranging from 41 up to 401. Data collapse, shown here 
for C(l), is obtained for v = 1, consistent with the model 
being in the universality class of the Ising model. In the inset 
is shown the divergence at the critical point for the infinite 
system. 

This figure contains the data for the lattice size rang- 
ing from N — 41 up to N = 2701. These results show 
that all the key ingredients of the finite size scaling are 
present in the concurrence. We note that finite size scal- 
ing scaling is fulfilled over a very broad range of values of 
N which are of interest in several protocols in quantum 
information. 

A similar analysis can be carried on for the next- 
nearest neighbour concurrence (7(2). Since d\C(2)\ x 
0. the logarithmic singularity here appears first in the sec- 
ond derivative with respect to A (see Fig.|| legend) . In the 
thermodynamic limit d\C(2) = 0.108 In |A — A c | + const. 
In this case also, the data collapse and finite size scal- 
ing (Fig.|J) agree with the expected scaling behaviour, 
v = 1. This completes the analysis of entanglement for 
the one-dimensional Ising model. 

A cornerstone in the theory of critical phenomena is 
the concept of universality - that is, the critical proper- 
ties depend only on the dimensionality of the system and 
the broken symmetry in the ordered phase. Universal- 
ity in the critical properties of entanglement was verified 
considering the properties of the family of models defined 
in Eq.(l) with 7 = 1/2. The range of entanglement £g is 
not universal. The maximum possible distance between 
entangled pairs increases and tends to infinity as 7 tends 
to zero. From the asymptotic behaviour of the reduced 
density matrix [|l9| we find that £g goes as 7 _1 . This 
however has no dramatic consequences; the total concur- 
rence C (n) stored in the chain is a weakly increasing 

function of 7 (for Q < 7 < 1, < En C ( n ) < °- 2 )- More 
interesting is the critical behaviour of the concurrence. 
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To be specific we consider C(l) in the case 7 = 0.5, 
shown in Fig. 4. As it was obtained for the Ising model, 
scaling is fulfilled with the critical exponent v = 1 in 
agreement with the universality hypothesis. 

The analysis of the 'resource' entanglement for a con- 
densed matter system close to a quantum critical point 
allows to characterize both quantitatively and qualita- 
tively the change in the wave-function of the ground state 
on passing the phase transition. A notable feature which 
emerges is that though the entanglement itself is not an 
indicator of the phase transition, an intimate connection 
exists between entanglement, scaling and universality. In 
a way this analysis allows us to discern what is genuinely 
quantum in a zero-temperature phase transition. The re- 
sults presented here might be tested by measuring differ- 
ent correlation functions, for example with neutron scat- 
tering, and extracting from these the entanglement prop- 
erties of the ground state close to the critical point Q . 

We finally discuss the perspectives of this work towards 
quantum computation. First, system sizes considered 
here (~ 10 3 ) could be those of a realistic quantum com- 
puter. Second, the scaling behaviour found could be a 
powerful tool to evaluate (and hence to use) entangle- 
ment in systems having different number of qubits. In 
particular, close to the critical point, the entanglement 
depends strongly on the field - so it could be tuned, re- 
alizing an 'entanglement switch'. Last, long range corre- 
lations, typical of the critical region, might be of great 
importance in stabilizing the system against errors due 
to imperfections. 
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